Goldstone modes and electromagnon fluctuations in the conical cycloid state of a multiferroic 
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Using a phenomenological Ginzburg-Landau theory for the magnetic conical cycloid state of a multiferroic, 
which has been recently reported in the cubic spinel CoC^CU, we discuss its low-energy fluctuation spectrum. 
We identify the Goldstone modes of the conical cycloidal order, and deduce their dispersion relations whose 
signature anisotropy in momentum space reflects the symmetries broken by the ordered state. We discuss the 
soft polarization fluctuations, the 'electromagnons', associated with these magnetic modes and make several 
experimental predictions which can be tested in neutron scattering and optical experiments. 
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I. INTRODUCTION 

Although feiTomagnetism and antiferromagnetism are the 
two most widely studied forms of magnetic order, more 
complicated, spatially modulated magnetic order parameters 
are also important and interesting from both fundamental 
and technological perspectives. A salient example, which 
occurs in the new class of ' multiferroic s ' ^i 3 ' 4 - materials 
that display an amazing coexistence and interplay of long 
range magnetic and ferroelectric orders - is magnetic trans- 
verse helical, or 'cycloidal', order. This order has acquired 
prominenco 3 - 4 - 5 ! 6 ! 7 ! 8 ! 9 ! 10 ! 11 ! 12 ! 13 ' 14 ! 15 ! 16 ! 17 ! 18 since it can induce, 
via broken spatial inversion symmetry 6,7 , a concomitant elec- 
tric polarization (P) in a class of ternary oxides, leading to 
interesting physics of competing and colluding ordering phe- 
nomena as well as potential applications 1 ' 2 ' 3 ' 4 . Among the ex- 
citing class of multiferroic materials, the cubic spinel oxide 
CoCr2C>4 is even more unusual, since it displays not only the 
coexistence of P with a spatially modulated magnetic order, 
but also with a uniform magnetization (M) 1 — in a so-called 
'conical cycloid' state (see below). 

Since in the conical cycloid state, the long range magnetic 
and polar orders are intertwined, it is crucial to understand the 
associated soft modes (i.e., low energy collective excitations), 
which should also be 'hybridized', leading to intriguing po- 
tential applications based on the electronic excitation of spin 
waves 19 and vice versa. A second motivation for studying the 
soft collective mode spectrum of a system with a complicated 
set of order parameters, such as the conical cycloid state, is 
that the Goldstone modes themselves caricature the underly- 
ing pattern of the broken symmetries, and thus, strengthen the 
understanding of the ordered state itself. In this paper, we do 
this by first identifying the magnetic Goldstone modes (i.e., 
magnons or spin waves) of the conical cycloidal order and 
deducing their dispersion relations which, as we clarify, sim- 
ply reflect the complex, anisotropic pattern of the underlying 
broken symmetries. We make several predictions for inelas- 
tic neutron scattering experiments based on our results for the 
magnetic fluctuations. We then identify the associated soft 
polarization fluctuations, which constitute a dielectric mani- 
festation of the magnetic modes, 'electromagnons', which can 



be observed in optical experiments. The interesting interplay 
of magnons and electromagnons in cubic multiferroics is the 
topic of this paper. 

CoCi"204, with the lattice structure of a cubic spinel, en- 
ters into a state with a uniform magnetization at a temperature 
T m = 93 K. Microscopically, the magnetization is of ferri- 
magnetic origin 16 , and in what follows we will only consider 
the ferromagnetic component, M, of the magnetization of a 
ferrimagnet. At a lower critical temperature, T c = 26 K, the 
system develops a spacial helical modulation of the magneti- 
zation in a plane transverse to the large uniform component. 
Such a state, for general helicoidal modulation transverse to 
the uniform magnetization, can be described by an order pa- 
rameter, 



M h = miei cos(q • r) + m 2 e 2 sin(q ■ r) + m 3 e 3 , (1) 



where {e^} form an orthonormal triad. When the pitch vector, 
q, is normal to the plane of the rotating components, the rotat- 
ing components form a conventional helix 2 -. A more compli- 
cated modulation arises when q lies in the plane of the rotat- 
ing components. For = 0, we will call such a state, which 
has been recently observed in a number of multiferroic ternary 
oxide o 3 ' 4 ' 7 ' 8 ' 9 ' 10 ' 11 ' 12 ' 13 ' 14 ' 15 , an 'ordinary cycloid' state because 
the profile of the magnetization resembles the shape of a cy- 
cloid. The cycloid state with ^ will be called a 'conical 
cycloid' state, because the tip of the magnetization falls on the 
edge of a cone, see Fig. Q] This is the low temperature mag- 
netic ground state in CoCr204, and is responsible for its many 
unusual properties, for e.g., the ability to tune P via tuning the 
uniform piece of the magnetization by a small magnetic field 
~ 0.5 T 16 . Notice that these states break the spin rotational 
and the coordinate space rotational, translational and inver- 
sion symmetries. It is easy to visualize that the helical, but 
not the cycloidal, modulation preserves a residual coordinate 
space U(l) symmetry (followed by a translation) about the 
pitch vector. 
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FIG. 1 : (Color online) The conical cycloid state with the mean field 
order parameter given in Eq.|2] The spins rotate along the pitch vec- 
tor, q, on the cycloidal (x — y) plane. The uniform component, 7713, 
which is along the z direction makes the magnetization tip fall on the 
edge of a cone. The polarization, P, is perpendicular to both q and 
the uniform magnetization. 



II. INTUITIVE UNDERSTANDING OF THE GOLDSTONE 
MODES 

To gain an intuitive understanding of the Goldstone modes, 
let's first consider the broken symmetries of the conical cy- 
cloid state, with a representative mean-field order parameter, 

M c (r) = (mi cos(gcc), ?7i2 sin((7a;), m 3 ), (2) 

shown in Fig. Q] As mentioned above, this state breaks the 
spin space rotation and the coordinate space rotation and 
translation symmetries. Note, however, that the translation 
symmetry is broken only in the direction of q. Since transla- 
tional symmetry is spontaneously broken in this system, uni- 
form translations along the direction of q, which can be pa- 
rameterized by the phase fluctuation ip(r), where the fluctu- 
ating magnetization may be given by M(r) = (mi cos(qx + 
<p(r)),TO2 sin(qx + ip(Y)),m, 3 ), must be a Goldstone mode. 
It is important to realize, however, that the elastic energy for 
this fluctuation cannot involve (d y ip) 2 , (d z (p) 2 , while it must 
involve the longitudinal component, (d x (p) 2 . This is because 
a uniform rotation of q, ip(r) = ay + (3z, rotating the pitch 
vector from (q, 0, 0) to (q, a, (J) must not cost any energy since 
the underlying Hamiltonian is assumed to be rotationally in- 
variant. The elastic energy must include (d x ip) 2 , however, 
since a change of the magnitude of q does cost energy. Thus, 
in momentum space, the dispersion relation for this Goldstone 
mode should be much softer in the directions transverse to q 
than in the longitudinal direction. 

The absence of a residual symmetry about q gives rise to 
a second Goldstone mode in the conical cycloid state. Notice 
that a uniform rotation of the cycloidal plane and the uniform 
magnetization about q does not cost energy, and therefore, 
such a rotation at long wavelengths must cost vanishing en- 
ergy. In the conventional helical state, this mode is already 
contained in the phase ip, since a uniform translation of a cir- 
cular helix long its pitch axis (i.e., a uniform (p) is equivalent 
to a rotation about the pitch axis by ip. The Goldstone mode 
fluctuations in the conical cycloid state are depicted pictorially 
in Fig. |2] 



III. GINZBURG-LANDAU HAMILTONIAN 

Since M and P respectively break time reversal and spa- 
tial inversion symmetry, the leading P-dependent piece in a 
Ginzburg-Landau Hamiltonian density, hp, for a centrosym- 
metric, time reversal invariant system with cubic symmetry 

is£, 

hp = P 2 /2 X + aP ■ M x V x M, (3) 

where \ > and a are coupling constants. We assume that 
P is a slave of M, in the sense that a non-zero P only occurs 
due to the spontaneous development of a magnetic state with 
a non-zero M x V x M. We consider a full Hamiltonian 
that is completely invariant under simultaneous rotations of 
positions and magnetization. This guarantees that any phase 
that can occur in our model is necessarily allowed in a crystal 
of any symmetry. The full Hamiltonian is given by 21 , H = 
J(h M + h P )dr = J hdr. Using P = -\aM x V x M 
to eliminate P, we can write the total Hamiltonian density h 
entirely in terms of M, 

h = tM 2 + 11M 4 + K (V • M) 2 + Kx (V X M) 2 
+K 2 M 2 (V ■ M) 2 + K 3 (M • V x M) 2 
+K 4 \M x V x M| 2 

+ J D L |V(V-M)| 2 + £ T |V(VxM)| 2 , (4) 

where we have u, Dl t > for stability. Due to competing 
magnetic interactions, some of the Ki can be negative. 

To discuss the parameter space for the conical cycloid state, 
t is assumed to cross zero at T m , and the system enters into 
a state with a uniform magnetization 7713 = \J—t/2u. As T 
drops further, the elliptic conical cycloid state, with the uni- 
form magnetization normal to the cycloidal plane and q in the 
plane of the cycloid, i.e., with a representative order parameter 
given by Eq. [2] is the lowest energy state in the regime t < 0, 
K 3 < 0, K < 0, < K 1 « -K 3 m\. In this regime, Eq.|3] 
defines the ground state among all the possible states with ar- 
bitrary mutual angles between the uniform magnetization, q, 
and the cycloid plane. K2, and K4 are relatively unimpor- 
tant for this state (Eq. |2] satisfies the saddle point equations 
with or without them), therefore, in what follows, we will set 
K2 = K4 = for simplicity 21 . 



IV. GOLDSTONE MODES IN THE CONICAL CYCLOID 
STATE 

To identify the Goldstone modes and to calculate their cor- 
relation functions, we follow standard methods: we first write 
M as its mean-field solution (describing the conical cycloid 
state) plus the fluctuations. We then substitute this total M 
in the Hamiltonian, Eq. |4] and expand the Hamiltonian to 
the second order in the fluctuation modes. A straightforward 
(though tedious) diagonalization of the fluctuation piece of 
the Hamiltonian would then produce the fluctuation modes 
(eigenvectors) and their energy dispersions (eigenvalues). As 



we will see below, there are four fluctuation mode 
ical cycloid state, among which two are massive ; 
two (a and S x , see below) are soft (Goldstone n 
long wavelength limit. By inverting the fiuctuatic 
Hamiltonian, one can also read-off the correlati 
of the soft modes from the matrix elements. 

Te begin, we write the total magnetization as 
5M, where <5M describes the fluctuations abov 
point solution M c . Generally, M can be written 

—m 3 8 y + mi cos (qx + </?)— m 2 S z sin 
M = | —m 3 S x + m 2 sin (qx + ip) + mi8 z cos 
7773 + Syirii cos (qx + ip) + S x m 2 sin i 

where ip describes the fluctuation of q, and S y am 
the rotation of the cycloidal plane and 777,3 a bo 
the z axes, respectively. Note that, for the circu 
state (mi = 7772), 5 Z can be taken to be zero 
renormalizes ip in this case. S x describes the ro 
cycloidal plane about the pitch vector itself. Exp 
first order in the fluctuation variables, we have 

-m 3 S y — (tpm\ + 6 Z 7772) sinqx 
SM = j —m 3 6 x + ((pm 2 + S z mi) cosqx 
S y mi cos qx + 8 x m,2 sin qx 



(6) 



To obtain the soft modes, we expand the Hamiltonian to 
second order in <5M. It is easy to check that the coefficient 
of the first order term is zero from the saddle point equations. 
The second order gives 



SH = 



t(6M.y - 

+D L |V(V -SM)\ 2 



2M 2 C (SM) 2 + 4 (M c • SM) 2 



•Z> r |V(V x SM) 



+K 



(V • 5MY 



(V x SM)' 



+K 3 [SM • V x M c + M c • V x SM) 2 

+2K 3 [M c • V x M c ] [SM ■ V x SM] (7) 

Substituting Eq. [6] into Eq. [7] taking the Fourier 
transform, and denoting Sq — ip, we find, SH = 
J2 P ,i,j P) r y(p)^(p)' where, i and j run from to 3. 
For brevity, we omit the full form of the 4x4 matrix V here. 

We should note, at this point, that in order for the fluctuation 
mode ip, and, in effect, the direction fluctuation of q to cost 
vanishing energy for infinite wavelengths, the cycloidal plane 
and the uniform magnetization themselves must rotate about 
the y and the z axes. The true Goldstone mode, for the third 
rotation fluctuation S x = 0, must then be a linear combination 
of ip, S y and S z . To capture this soft mode, we first take S x = 
and diagonalize the resulting 3x3 matrix. The eigenvalues 
for two eigenvectors, (3, 7, remain non-zero even when the 
momentum p — » (massive modes), but the other eigenvalue 
becomes zero in this limit (soft mode). The corresponding 
eigenstate of the soft mode, to linear order in p — |p|, is given 

by, 





FIG. 2: (Color online) Left: The magnetic Goldstone mode a of the 
conical cycloid state. For an arbitrary small fluctuation of the pitch 
vector, q, the cycloidal plane and the uniform magnetization must 
also rotate by angles S y and S z about the axes y and z, respectively, 
for a zero energy deformation. Right: The other Goldstone mode S x . 
S x describes the rotation fluctuation of the entire system about the 
pitch vector, and, if spatially uniform, costs no energy. For clarity, 
the rotations of 7713 and P are not shown. 



This is one of the two cycloidal Goldstone modes found in 
this paper, see Fig. [2] To order p 2 , we have the corresponding 
eigenvalue, 



= 2(7 



2 2 

q p x - 



(9) 



As expected, there is no contribution from p y ,p z at this or- 
der. As emphasized before, this is a reflection of the rotational 
symmetry of the underlying Hamiltonian. The next higher or- 
der contribution to the Goldstone mode eigenvalue is given by 
Wi = uip A y + u 2 pi + u 3 p x + u 4 plpl + u b p 2 y p 2 z + u 6 plp 2 z , 
where the m's are functions of 777,1, "72, "73 and the coupling 
constants Kq, K\, K 3 , Dl t- 

The other Goldstone mode of the conical cycloid state is 
simply the mode S x , see Fig. [2 with the momentum space dis- 
persion relation starting at the order p 2 ,Py,p^- As explained 
before, spatially uniform rotation of the whole system about 
the direction q = x does not cost energy, so the long wave- 
length fluctuations, represented by S x , cost vanishing energy. 

In the presence of lattice and spin anisotropics, the forego- 
ing results are valid only above the anisotropy energies. The 
anisotropic dispersion of the mode a crosses over to a more 
isotropic dispersion, one which depends quadratic ally on all 
of p x ,p y ,p z , below the lattice anisotropy energy. However, 
it continues to remain a true Goldstone mode because of the 
broken translational symmetry. In this respect, this cycloidal 
magnon is analogous to the phonon mode in a crystal, rather 
than a true magnon mode. Below the weak spin anisotropy 
energy, the other Goldstone mode, S x , should acquire a gap 
given by this spin anisotropy energy. 

In the most general case, the two soft modes will couple. In 
terms of the corresponding eigenstates, the 4x4 matrix can 
be rewritten as a 2 x 2 matrix (plus unimportant contributions 
coming from the massive modes), 

2 (777^ + 777^) q 2 p 2 x + ui -PxPzVo ~ iPxPyPzVl/q 
-PxPzVo + ipxPyPzVi/q f{p 2 ) + m|D T p 4 /2 

(10) 

where vo, v\ are constants and f(p 2 ) is a second order poly- 
nomial function of p x ,p y ,p z . By inverting this matrix, we 
find, 



ip (p) + ip z S y (p) /q + ip y S z (p) /q. 



(8) 



C aa (p) = V ffcl>?/A(p) 

' * ■} . — T II K 



t—x,y,z 
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Cs^sAp) = [4 (.9i + 92) q 2 pl + wi] /A (p) 

Ca<5x (P) = PxP Z «o/A(p) (11) 

where C M „ (p) = {fi (-p) f (p)), A (p) 

PlT,i= x ,y,zl 3 iPi + w iE != x, 9 ,zW? + -■ is the deter- 
minant of the matrix ( TfOl i. and the rji 's and the /Jj's are 

constants. Remarkably, forp^ = 0, we find, 

C aSx (p) = 0, (12) 

so there is no contribution from p y and p z to order p 2 in the 
C aa (p) correlator, as expected. 



V. MAGNETIZATION CORRELATIONS AND NEUTRON 
SCATTERING 

From the energy resolved neutron scattering cross sections 
near p = q, it should be possible to track the p-space disper- 
sions of the fluctuation modes a, (3, 7 and 8j- A . Most notably, 
the anisotropic dispersion ~ uj + u>i of the mode a, carica- 
turing the complex broken symmetries of the conical cycloidal 
order, should be experimentally testable. 

Using the soft mode eigenvectors, we can calculate 
the full static magnetic susceptibility tensor, Xij — 
{Mi (-p) Mj (p)) = (SMi (-p) SMj (p)). For instance, the 
dominant terms of xu are 

p 2 1 

Xxx ~ "T,3^C QQ (p) + -TO 2 (C Qa (p-q)+C Qa (p + q)) 
q z 4 

Xyy ~ m\Cs m s x (p) + -fnl (C aa (p - q) + C aa (p + q)) 

1 1 

Xzz ~ ^5 £ (C Q a(p-q) + C QQ (p + q)) 

+ \™l (C SmSm (p - q) + Cs.6. (P + q)) 

-P±m im2 ( f" s * + q ) - ^ ( p - q \ ) (13) 
4q V +C*.o(p + q) -Ca ia (p-q) 7 V 

It follows that the susceptibility functions diverge both at p = 
and p = ±q for the conical cycloid state, the divergence at 
p = originating from the fluctuations of ■ 

The susceptibility functions show different behaviors when 
p approaches ±q or along different directions in momentum 
space. For instance, when p — » q along p x , all xu diverge as 
(Px — q)~ 2 - On the other hand, when p — > q along y or z 
directions, Xxx and Xyy scale as p~ A , and Xzz scales as p~ 2 . 
In neutron scattering experiments, the following quantity is 
related to the frequency integrated scattering cross section 22 , 

f°° , 1 / / u)\\ d 2 a 



where x (p) — y>ij ~ ^p"J Xij> u is the frequency and fi is 
a solid angle. Near p = q, the dominant terms in x (p) are, 

X (P) ~ jrnlC aa (p - q) + \m\C &x&x (p - q) . (15) 

When p — > q along p x , x (p) ~ (p x — l)~ ■ In contrast, 
when p — > q from the p y or p z directions, the divergence 
goes as x (p) ~ Pi*- 

VI. POLARIZATION CORRELATIONS AND 
ELECTROMAGNONS 

The static dielectric susceptibility tensor, Xijt is propor- 
tional to the polarization correlation functions, Xy(p) x 
(Pi(— p)Pj(p)). They can be straightforwardly derived by 
using P~MxVxM and the magnon correlation functions 
C^v (p). For brevity, we do not give here the full expressions 
for the polarization correlation functions. Typically, the corre- 
lation functions transverse to P diverge near p = and p = q 
due to the magnetic Goldstone modes in the conical cycloid 
state. Since the underlying magnons manifest themselves in 
the dielectric response of the system, these fluctuations are 
sometimes called 'electromagnon' fluctuations^ 2 ^. 

Since the typical optical wavelengths ~ 0(100 nm) are 
much longer than the lattice constants ~ O(IA), we only dis- 
cuss here the behavior near p ~ 0. Note that the fluctuations 
near q may also be influenced by the so-called symmetric cou- 
plings between P and M 24 , which do not contribute to the 
uniform macroscopic P. We will ignore these effects here 
since they are not accessible by the experiments. The trans- 
verse correlator along the direction of 7713 always diverges 
in this limit, (P z (-p) P z (p)) ~ pr 2 - (i = x,y,z,pj& 
Pi) = 0). This divergence arises from the mode S x , which 
rotates the cycloidal plane about x yielding a fluctuation of 
P along z. The other transverse susceptibility also diverges, 
(P x (-p) P x (p)) ~ p^ 2 , for p x ,p z = 0. This divergence 
arises from the Goldstone mode a. Note that the mode a 
includes the rotation fluctuation S z , which induces a polar- 
ization fluctuation along x. These characteristic divergences 
should be observable as peaks in the appropriate static dielec- 
tric constants, revealing the existence of the electromagnon 
fluctuations in the conical cycloid state. In the conical cycloid 
state, but not in the ordinary cycloid state, the polarization 
correlation functions diverge also near p = q, the coefficient 
of proportionality of the diverging piece being 7713, but these 
electromagnon fluctuations will be difficult to see in optical 
experiments because of the non-zero momentum. 

VII. CONCLUSION 

To summarize, we have identified and discussed the mag- 
netic and polarization fluctuation modes of the conical cy- 
cloidal order in a multiferroic. One of our primary predictions 
is the unusual dispersion relations of these soft modes, which 
can be experimentally tested on CoCr204, thereby revealing 
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the complex pattern of the broken symmetries and their asso- 
ciated Goldstone modes. We also predict the divergence of the 
magnetization and the polarization correlation functions; the 
latter reveals the hybridized soft mode, the electromagnon. 
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